Effect of Short-Range Fluctuations on Thermodynamic and Resistive Properties: The 

case of Ising Order 

Mukul S. Laatfl 

Institut fur Theoretische Physik, RWTH Aachen University, 52056 Aachen, Germany 



o 

> 

O 

O 

li 
o 
o 

u 

Q-i! 



c 

o 
o 



> 

00 
On 
O 

d 
o 



X 



Luis Craco 

Technische Universitat Dresden, Institut fur Physikalische Chemie und Elektrochemie, 01062 Dresden, Germany 

(Dated: October 7, 2010) 

We consider the effects of the non-local Ising-like "core spin" correlations on the order-parameter 
fluctuation contribution to the resistivity and thermodynamics of metals showing Ising-like order 
at finite temperature. We employ the well-known cluster-variation method, and present explicit 
results in the pair approximation for short-range order. Our calculation generalizes earlier works, 
where such effects were considered in the mean-field (Ornstein-Zernicke) approximation. The mean- 
field (MF) transition temperature T^ rF , is corrected to 0(1/(2), and the effect of the Ising spin 
fluctuations on the dc resistivity and magnetothermal responses is analyzed in detail. The method 
can be applied straightforwardly to lattices in arbitrary d, and, as an appealing feature, it reproduces 
the exact correlation length and T c ld = for the Id Ising model. We apply our results for two 
interesting physical cases: (i) the double-exchange model with Jh » t, where the core-spins can 
be approximated quite well by Ising spins, and (ii) a model of band electrons coupled to a localized 
subsystem which undergoes a nematic ordering transition coupled to an appropriate structural 
transition. 

PACS numbers: 71.28+d,71.30+h,72.10-d 



I. INTRODUCTION 

Finite temperature (T) order in solid state systems 
leave their imprint on thermodynamic and transport re- 
sponses across the order-disorder transition. The steady 
enhancement of the correlation length of the system as 
T is lowered toward T c from above directly impinges, in 
particular, on the T-dependent carrier scattering rate, 
and manifests itself as changes in thermal and transport 
behavior as the phase transition is approached. For a 
classical transition from a ferromagnetic metal to a para- 
magnet, such effects on the dc resistivity for the Heisen- 
berg universality class were considered more than forty 
years ago^. This estimation has been invaluable in inter- 
preting a large body of data on band ferromagnets. 

Many real compounds, however, show pronounced ef- 
fects of short-range order (SRO) on carrier properties. 
Devising computational schemes capable of incorporat- 
ing such SRO effects into realistic microscopic calcula- 
tions has long been the goal in the field of electronic 
structure. However, to date, it has proved to be ex- 
tremely demanding to go beyond the single-site approxi- 
mation: in disordered transition-metal alloys, the coher- 
ent potential approximation (CPA) has long been used 
with great success^, while for strongly correlated elec- 
tronic systems, dynamical-mean-field-theory (DMFT) in 
conjunction with ab initio (LDA) band structure also 
holds a lot of promise 3 . However, by construction, such 
a single-site approach cannot capture the full effects of 
SRO on carrier scattering rates. This requires extending 
DMFT(CPA) to capture SRO effects via proper cluster 
extensions, a formidable task realizable currently only 
for one-orbital models 3 . Moreover, details of such SRO 



induced effects on carrier scattering rates and transport 
properties remain unexplored avenues for study, to our 
best knowledge. 



Motivated by the above, we focus on a much "simpler" 
version of the above problem that nevertheless bears on a 
number of interesting issues in modern condensed matter 
physics. Specifically, we consider physical systems with 
an Ising like phase transition to an ordered phase (either 
rigorously or to a very good approximation). It turns 
out that, for this particular case, a controlled estimation 
of corresponding order parameter fluctuations to 0(1/ d) 
(with d the system dimensionality) and their feedback 
on band-like carriers can indeed be carried out under 
restricted conditions (see below). 



To emphasize the relevance of our "simple" theory for 
physical problems of great current interest, we consider 
two explicit cases: (i) the double exchange (DE) model 
with Jh >> t, in which case, the spin correlations can 
be well approximated by an effective Ising model. This 
is believed to be an appropriate low-energy model for 
certain colossal magnetoresistance manganites in a lim- 
ited doping region corresponding to a well-ordered low 
T ferromagnets. (ii) the problem of an electronic system 
undergoing an electronic nematic transition^, coupled, by 
symmetry, to an appropriate lattice distortion. This or- 
dering rigorously belongs to the universality class of an 
Ising model in an external zeeman field. Evidence for 
electronic nematicity appears in various cases of great 
current interest, like bilayer SrzR^O-j, underdoped Fe 
arsenides and underdoped high-T c cuprates. 
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II. FORMALISM 

A. Boundaries of Validity 

In this section, we describe the well-known cluster 
variation method (CVM) for the simple DE model ((i) 
above) . But obviously, the formalism applies equally well 
for case (ii) with replacement of the real Ising spin by a 
5=1/2 pseudospin corresponding to Ising nematic or- 
der. 

We start with the double exchange (DE) model in the 
limit >> t with a nearest neighbor ferromagnetic 
(FM) coupling between core spins^. This model is be- 
lieved to be applicable to doped Europium oxide (EuO) 
and is an effective low-energy model for well-doped man- 
ganites. It reads 

h = J2 MS)(c 4 W+M-J' s?s?-h ext J2si • 

<ij><? <ij> i 

(1) 

In the above, we have included only the Ising part of 
the intersite FM coupling, because the transverse spin 
fluctuations make a negligible contribution when Jh » 
t. In this limit, the spin correlations are well approxi- 
mated by that of an effective Ising model, with an effec- 
tive intersite coupling J' = (t 2 / Jr)x{1 — x) for the well- 
doped FM metal. Also, %(S) = t[l+ < S Z S Z > /2S* 2 ] 1/2 
and h ex t is the external magnetic field. We emphasize 
that tij is to be understood as a renormalized carrier 
hopping integral, reduced from its bare band structure 
value by strong local Hubbard correlations in the real 
system. 

In the case of the nematic order, the "core spin" above 
is replaced by a nematic pseudospin, N z , which quantifies 
the degree of electronic anisotropy, reflecting spontaneous 
breaking of (discrete) lattice rotational symmetry 5 . In an 
Ising nematic, this is defined simply as N z = 
where n a is the average fermion occupation number along 
a(= x,y). Physically, this can arise from multi-orbital 
effects (as probably the case in Sr 3 Ru20i^ and under- 
doped Fe arsenides^) or from anisotropic (effective) two- 
body interactions^ (as possibly in underdoped cuprates). 
We will assume that electronic nematic (e-nematic) order 
has occured, without specifying its microscopic origin in 
this work. By symmetry, N z directly couples to the lat- 
tice strain, which acts as a field term conjugate to the 
order parameter. The Hamiltonian is 

H= £ t ij (clc ja +h.c)-J'Y / NfN?-5(T)Y,Nf , 

<ij>cr <ij> i 

(2) 

where 8(T) is the lattice strain, simplified here to a single 
T-dependent function to illustrate its dominant role. The 
intersite exchange is "ferro", since the e-nematic state 
involves a q — instability. Again, tij is to be understood 
as an effective hopping of (substantially) renormalised 
band carriers in reality. 



Effects of arbitrary SRO on carrier dynamics in real 
materials has hitherto relied on heavy-duty numerical ap- 
proaches*^. Ideally, one would like to extend the local 
approximation in a way that is amenable to treatment 
of the problem in any dimension, and is computationally 
tractable. In what follows, we present an analytical ver- 
sion of the cluster- variation method (CVM), which has 
been very widely used in studies of alloy formation and 
stability. This scheme satisfies the above requirements, is 
physically transparent and amenable to easy numerical 
implementation. We emphasize that, as a first step, we 
look at these effects within effective models, and our ap- 
proach is only valid in the cases where SRO arises from 
(electronic) degrees of freedom that are effectively decou- 
pled from the band-like carriers at low energy. This is 
also the reason why we cannot describe the microscopies 
of the ordering mechanisms: we take this as given. 



B. Cluster Variation Method 

As a first step, we realize that the core (Ising) spins can 
be thought of as an alloy system of 5 2 = —1/2, +1/2. A 
perfect FM order (low T) corresponds to a pure system 
of either S z = -1/2 or +1/2. At sufficiently high T, the 
system is non-magnetic and corresponds to a completely 
random alloy of S z = -1/2, +1/2, where the CPA should 
be a very good approximation. At intermediate T, how- 
ever, SRO effects are dominant. In this situation, we 
map our problem to that of a short-range ordered binary 
alloy, which has been extensively studied in the field of 
magnetism of disordered transition metals. In particular, 
we employ the CVMii, which has been very successfully 
used in this context. The CVM gives analytical expres- 
sions for the correlation functions of the system within 
the so-called pair and square approximations^ 2 -. In addi- 
tion, it satisfies the "diffuse intensity sum rule" exactly 
in d — 1, a feature that is not shared by other, less so- 
phisticated approximations. We make the identification 
Sf = -1/2 ->■ n, = and S z = 1/2 -> n, = 1. We start 
by writing down the functional: 

F[a] = E(a)X(a) + k B Tj2x(a)lnX(a) , (3) 

a a 

where (for a binary system such as the one we consider) 
for N lattice sites, the sum runs over 2 N configurations. 
Here, E(a) and X{a) denote the energy and probability 
of a configuration a. The system free energy is just the 
minimum of F[a], where the minimization is carried out 
over all density matrices X(a) subject to normalization: 

In general, the average configurational energy of the 
spin subsystem can be written down as a linear combi- 
nation of many-body correlations. For a binary system, 

E{a) = 1/2 J'jSlS z + £ hiS? , (4) 

ij i 
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and 



F[a] = E(a) + -^— y<T<j'{h3)lny„v'{i,j) 



k B T(l - z)) j x a (i)lnx<j{i) , 



(5) 



where x a (i) = [1 + cr£i(«)] is the single-site occupation 
probability and yaa'{hj) = (1/4) [1 + cr£i(i) + + 
j) + (ycr'&ihj) is the pair probability. The quantities & 
with i = 1,2 are given by derivatives of the free energy 
w.r.t the inhomogeneous field: =< Sf >= 4jE- and 



<*ij = Uhj) =< s*s* > - < Sf >< s? >= 



d 2 F 
dhidhj 



(6) 

and the spin susceptibility including SRO contributions 
is calculated as the Fourier transform of . 

In what follows, we will consider only the pair approxi- 
mation^ as a first step, so that (Ri — Rj) = a, the lattice 
constant. To compute the thermodynamic response and 
carrier scattering rate, we need the free energy of an ef- 
fective Ising model in a zeeman field: interestingly, in the 
pair approximation, we show that these can be obtained 
analytically. We find the free energy of the "core-spin" 
subsystem in the pair approximation as 



F(m, a) = N(-zJ'— - mh ext ) - zJ'(l - m 2 )a + F[S] , 



where 



F[S] 



(7) 



k B T 



ksTz 



m 1 
— In — 



m 1 
— In — 



a 



(8) 



The magnetization, m(T), is computed as the first 
field-derivative of the free energy: 



m(T) = tanh 



zJ'm(l — a) + h 



k B T 



(9) 



where a = a i3 = [< S!f5? > -m 2 (T)]. Notice the differ- 
ence from the Weiss mean-field (MF) result; the non-local 
spin correlation function explicitly enters the transcen- 
dental equation for m(T). 

The q-dependent core-spin correlation function is de- 
rived from the second field-derivative of the free-energy: 

x(q> T) = ^^ e ^(R.-Rj)[< S?S* > - m 2 (T)} (10) 

id 



X 1 (q, T) — for q = 0; i.e, as the solution of 
l-(z- l)tanh(J' (3 C ) = 0. In d = 1, this gives T C FM = 0, 
and the correlation length, £(T) = tanh(J' f3), which 
coincides with the result known from the exact solu- 
tion^. It is easy to see that for small q around q = 0, 
X(q, T) ~ A{T)/^- 2 + q 2 x + ql + ql) in d = 3, which is the 
classical Ornstein-Zernicke behavior with the FM spin 
correlation length given by £, 2 {T) = 2tanh(J'/3)/[l — (z — 
l)tanh 2 (J' 0)\, again noticeably different from the MF 
result. Comparing with the MF result, we see that in- 
clusion of non-local FM spin correlations within the pair- 
approximation (CVM) depresses the transition tempera- 
ture: T C FM = T™ F /(z - 1), where T™ F = 2J'z/k B is 
the Weiss MF result, as expected. The q = susceptibil- 
ity is directly written as, 



X%T) = {g^ B ) 1 - (z l)tanh(J'P) 



(12) 



which deviates noticeably from the MF-Curie- Weiss form 
at intermediate T. 



III. SOME PHYSICAL OBSERVABLES: 
RESISTIVITY AND THERMODYNAMICS 



Given the spin correlation function, the spin disor- 
der contribution to the dc resistivity is evaluated using 
Fermi's golden rule from the equation^, extended to an 
Ising transition. 

Pdc(T) = — / X (q, T)(l - cos9)de (13) 
ne z J FBZ 

for D dimensions. Here, sin(9/2) = q/2kp, whith 
t _1 (T) = J FBZ x( ( i,T)(l — cos6)d6 again involving the 
non-local spin correlations via x(q, T). We expect the 
detailed nature of the field induced changes in Pdc(T) to 
be determined by the field-induced changes in the ther- 
mal order parameter fluctuations. An external field will 
increase carrier mobility via tij, reducing r~ 1 (T), and 
resulting in increased itinerance above T c . 

Calculation of the thermal expansion and magnetovol- 
ume effects requires a bit more work. In a system with 
Ising-like local moments (magnetic or nematic), which 
are stable as T is raised, one invokes the same local- 
ized picture as for insulators, where any magnetovolume 
term is due to the volume dependence of inter-site in- 
teractions. Following^, the order parameter fluctuation 
contribution to the thermal expansion coefficient is, 



*(q,T) = 



[l-Cl(T)] 



l + (*-l)£ 2 2 (T)-26(T) e ? 



(11) 



in d dimensions and = Y^a=i C0S [la{Ria — Rja)]- 
In general, T C FM is determined from the solution of 



(T\ K d 



ij 



Jl'zC Sj ^> 



! (T)] 



(14) 

where 7 m = —dlnJ/dlnV is the magnetic or lattice 
Griineisen parameter. The integrated magnetovolume is 
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then, 

T i =f^r <iT ^(5 j[<s ' s j > -" i2(T)i ) ■ 

(15) 

Using Maxwell's thermodynamic relation, the entropy 
change in a field can be directly computed from, 

AS M (T, h) = S M (T, h) - S M (T, 0) = / (—) h dh , 

Jo dT 

(16) 

using the equation for the order parameter in the 
pair approximation derived above. (A large change in 
ASm(T,H) is expected near Tf M in DE ferromagnets, 
as also in disordered local moment ferromagnets^ and 
would be interesting in the context of applications to 
magnetic refrigeration) . 

The above equations show how short-range order di- 
rectly influences various transport and thermodynamic 
quantities in metallic systems showing an Ising-like semi- 
classical order via coupling of band carriers to order pa- 
rameter thermal fluctuations. Obviously, our scheme is 
only valid when the Ising-like order arises from micro- 
scopic (electronic) degrees of freedom that are effectively 
decoupled from the band-like carriers, but, once estab- 
lished, affects carrier dynamics as shown above. 



IV. RESULTS 




FIG. 1: Order parameter as a function of T. The MF T c is 
reduced by 30 percent by order parameter fluctuations. In a 
finite "field" , clear hysteresis (Ising model in a zeeman field) is 
seen for moderate h. In the DE model, a small magnetic field 
results in appreciable FM short-range order above T c , modify- 
ing thermodynamic and resistive behaviors (see Figs l3l5 |). In 
the nematic case, as in underdoped Fe arsenides, the "field" 
is the strain,which acts as a conjugate field to the nematicity, 
i.e, it couples linearly to the nematic pseudospin. Nematic 
correlations persist well above Tt-o, in qualitative accord 
with recent indications in underdoped 122-Fe arsenides. 



We now discuss how the above "simple" formalism cap- 
tures, surprisingly, a wide range of interesting features in 
the cases (i) and (ii) mentioned in the Introduction. 



A. Colossal Magnetoresistance Manganites 

We note that the Langer-Fisher formulation has been 
applied with good success by Majumdar and Little- 
wooded. They considered a variety of doped magnetic 
semiconductors. However, effects of magnetic short- 
range order have received scantier attention in this con- 
text, and this is the issue we will address here. 

In this section, we show the results of our calculation 
and discuss them in some detail. In Fig. [TJ we show the 
core-spin magnetization as a function of temperature, T. 
Notice the difference from the MF result, where m(T) 
vanishes at Tf^ F (the MF curve is shown by the dot- 
ted curve in Fig. [T]). Magnetic SRO (J' > 0) smears 
the MF transition, as shown in Fig. [TJ the tail showing 
persistent FM SRO above Tf M . A small magnetic field, 
hext = 0.05t, aligns these short-range-ordered regions, 
giving an appreciable M(T) well above Tjf M . Obviously, 
pre-existing SRO, in conjunction with a small J' ~ t 2 / Jh 
makes it easier for a small h ext to polarize the high-T 
phase. Similar behavior is seen for h ext = O.lt. We 



also clearly see emergence of hysteresis in m(T), a fea- 
ture expected in an Ising model in a zeeman field. The 
static (q = 0) spin susceptibility, x(®>T) shows a dis- 
tinctly non-Curie- Weiss behavior below T ~ 2Tf M , as 
indeed observed experimentally, showing the persistence 
of magnetic SRO to rather high temperatures. 

The dc resistivity, pdc{T), computed within the Born 
approximation, actually exhibits an insulator-like behav- 
ior above T C FM (Fig. g]). Within simple MF theory, 
PdciT) = const above Tf M (see Fig.[2]in the case of clas- 
sical (or Ising) spins, since these act like non-magnetic 
scatterers, giving a resistivity characteristic of impurity 
scattering. With J' = 0, we thus obtain a result resem- 
bling that obtained from MF theory, as seen in Fig. [5] 
With J' ^ 0, this behavior is drastically modified: a 
sharp peak, reflecting coupling of carriers to the increas- 
ingly singular x(q, T) , is clearly seen in the inset of Fig. [1] 
We thus conclude that it is precisely the intersite corre- 
lations, and, in particular, the dominant FM SRO above 
Tf M which drives the insulator-like resistivity in the PM 
state. The effect of the external magnetic field is striking: 
hext — 0-05t drives the high T resistivity metallic. This 
is directly related to the field-dependent suppression of 
the FM short-range correlations, as seen in Fig. [3l re- 
ducing the scattering rate, enhancing the carrier mobil- 
ity and driving the system metallic. Similar (enhanced) 
metallic behavior is seen for h ex t = O.lt, and the dc re- 
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T 



FIG. 2: dc resistivity with J' = 0. Notice how the resistiv- 
ity shows a broad and smmoth change across T c . This is the 
expected result for the Ising case, which corresponds to "po- 
tential" non-magnetic disorder, giving Pdc(T) ~ const above 
T c . This is an artifact of the MFA and is corrected by the 
CVM (See Fig.0. 



sistivity shows a huge drop in pdc{T) as a function of 
h e xt- We emphasize that while a large MR drop in h ex t 
is plausible within MF treatments of the DE model with 
additional randomness (disordered static JT distortions, 
static chemical disorder), such scenarios require a large 
disorder (sufficient to split the doping-induced impurity 
band from the lower exchange-split band), the origin of 
which is unclear. Our study should describe Tl^M n2<97, 
where the local moments and band electrons arise from 
distinct orbital states, and do not hybridize much with 
each other—. 

We consider next the magnetic entropy change, com- 
puted from the T-derivative of the magnetization. In 
Fig. H we show S M {T,h ext ) for h ext = 0,0.05,0. It. A 
number of interesting features are apparent; the field- 
induced redistribution of Sm is drastic, the sharp peak 
around T C FM changes to a broad peak like feature in 
an external field. The quantity A 5a/ = [5a/ (T,H) — 
Sm(T,0)]/Sm(T,0), which measures field induced en- 
tropy change shows very encouraging behavior, attaining 
values up to 6 in h ext — 0.15i. We also notice that the 
overall shape of the curve (shown in the inset of Fig. [3} 
is in good qualitative agreement with observations^. In 
this context, we notice that Gd metal, widely used as a 
magnetic refrigerant, shows A5a/ = 4.0 near room tem- 
perature. Thus, Jh » t (small J) is one of the con- 
ditions favoring a large A5a/, and possible application 
of such systems to magnetic refrigeration. Further, the 
manganites exhibit considerably small magnetic hystere- 
sis with a coercivity of 50 Oe near T C FM , which should 
enhance their magnetic cooling efficiency. Finally, inter- 
estingly, we remark that the A5a/ (T) we extract qualita- 



18.0 




0.0 0.2 0.4 0.6 

T 



FIG. 3: dc resistivity as a function of T, with h = (black 
curve) and in a finite "magnetic" field (red and blue curves). 
Clear hysteretic behavior for intermediate h is obvious, and 
is related to that in m(T, h) vs h in Fig. [T] Clear qualita- 
tive differences from the results based upon Weiss mean-field 
approximation (MFA) (see Fig. [5J are manifest, and are con- 
sequences of the interplay between appreciable short-range 
order and its suppression by an external field. 

tively resembles that found for ternary CoMnGei-^Sna; 
ferromagnetic alloys by Hamer et al^. Our formalism, 
based on magnetic SRO, can also be readily extended 
to include structural/chemical SRO in disordered ferro- 
magnetic alloys. This entails incorporation of specific 
electronic structure details, and lies out of scope of the 
present phenomenological theory. 

Finally, we consider the magnetic contribution to the 
volume expansion, AVm{T)/V. Before presenting our 
results, we present a brief physical picture of magnetovol- 
ume effects in magnetic systems. At high-T, with com- 
pletely disordered local moments, one has as many t spins 
as 4,-spins, so the magnetovolume interactions average to 
zero. As T decreases, short-range local moment correla- 
tions begin to develop, and the magnetovolume interac- 
tion is determined by the spatial dependence of the spin 
correlation function. Obviously, the MF approximation 
would give incorrect estimates of magnetovolume effects 
in manganites, since, as we have seen, the FM spin cor- 
relations are very different from those expected from MF 
approaches. In view of the ability of the CVM to yield 
a more consistent description of these correlations, one 
expects that it is able to better describe magnetovolume 
effects in such systems. 

In Fig. [SJ we show the magnetovolume changes for dif- 
ferent field values. Note that it tends to zero for suf- 
ficiently high-T, consistent with the picture above. In- 
terestingly, at temperatures much higher than Tf M , the 
effect of short-ranged FM spin correlations is manifested 
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FIG. 4: Order parameter fluctuation contribution to the 
entropy as a function of T. Significant improvement over 
the mean-field estimate is clear. The field-induced entropy 
change, ASm(T), attains its maximum in a range of field val- 
ues where m(T, h) vs h shows hysteretic behavior. For "local 
moment" metallic ferromagnets, combination with a coerciv- 
ity of O(50) Oe would make them attractive candidates for 
magnetic cooling applications. Similar features have recently 
been reported across the T-O structural transition in under- 
doped Fe arsenides. 



FIG. 5: Volume change as a function of T. There is a sign 
change at the ordering temperature (T c ). Notice the negative 
sign of (AV(T)/V) below T c , followed by a non- linear vari- 
ation and hysteresis above T c when J' > 0. Similar features 
are characteristic of some Invar materials. They have also 
been reported recently across the T-O structural transition in 
underdoped Fe arsenides. 



B. e-NEMATIC ORDER IN METALS 



in cxm (T), which is negative and decreases with decreas- 
ing T upto about Tjf M . It is interesting to observe that 
a negative AVm(T)/V (thermal expansion coefficient) is 
characteristic of invar alloys^, and is observed both in 
the chemically disordered and ordered cases, demonstrat- 
ing that chemical disorder is not the driving force of invar 
effects. But thermal spin fluctuations introduce spin dis- 
order which could produce such effects. In our case, the 
magnetovolume changes sign at Tf , so more work is 
needed to describe such systems, which may include cases 
where there is no correlation between Tf M and the sign 
change in q>m(T). 

The field-dependence of AVm(T)/V shows how it is 
sensitively affected by h ext . The fractional change in 
AVm(T, h ex t)/V is large, and can decrease by as much 
as 175 percent in h ext = O.lt. More interestingly, it can 
be switched in sign by an applied field above Tf M , as is 
clearly seen from Fig. [5j This fact might be potentially 
useful for technological applications. The change is max- 
imum just around Tf M , and is thus intimately related to 
the field-induced suppression of spin fluctuations above 
Tf AI . Hence, it can only be accessed in a theory which 
describes SRO in a consistent way beyond the single-site 
limit. 



Here, we will show that, remarkably enough, our "sim- 
ple" technique also provides a simple qualitative explana- 
tion for a variety of ill-understood fluctuation effects in 
physical cases where an electronic nematic (e-nematic) 
instability arises in metals at low T. The whole struc- 
ture of the theory remains intact upon re-identifying the 
S z as nematic pseudospin variables, N z . Moreover the 
electronic nematic transition rigorously falls into the uni- 
versality class of an Ising ferromagnet in a zeeman field . 
Obviously, we want to emphasize that, in our effective 
approach, deeper questions relating to microscopic origin 
of such nematic state(s) cannot be answered: our main 
focus is to try to understand its consequences for thermo- 
dynamic and transport responses. To this end, we will 
also drastically simplify the more complicated situation 
in real systems of interest (like SrsR^O^ and under- 
doped Fe arsenides^—) by replacing their undoubtedly 
strongly correlated (multi) bands by a single free-electron 
like band. In such systems, the N z may microscopically 
arise via orbital selective Mott localization of a subset of 
the d orbital band states in the real correlated system, 
and, in our effective model, we assume that this has al- 
ready occured. Thus, our analysis cannot hold for study- 
ing nematic correlations in underdoped cuprates, where 
the possible nematic order and the "itinerant" carriers 
arise from a same single ban d 22 ' 23 . 

(i) First, 1/d (classical) order parameter fluctuations 
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depress T nem by about 30 percent compared to the Weiss 
MF result. Without coupling to strain, the transition is 
continuous, but turns into a first-order one with a finite 
(not too large) h ext — S(T), with hystcretic behavior. 
This is clearly seen in Fig. [3j where m(T) is now the 
nematic order parameter. The lattice constants will di- 
rectly react to the onset of nematic order, inducing a 
structural transition (tetragonal-orthorhombic in the Fe 
arsenides). With reference to underdoped Fe arsenides, 
this describes the tetragonal-orthorhombic (T-O) phase 
transition: this can be either a first- or second order tran- 
sition as we traverse the families of real Fe arsenides^. 

(ii) Correspondingly, the dc resistivity shows a clear, 
sharp anomaly at T nem for 8 = 0, and a "rounding off" 
of this sharp feature for a finite S, with a hysteresis in 
Pdc(T) above T nem , as shown in Fig. [3] 

(iii) The entropy change, AS nem , now interpreted as 
the entropy loss across the e-nematic transition, also ex- 
hibits a clear, relatively sharp, maximum precisely at 
Tnem-, as shown in Fig. 01 

(iv) Finally, the nematic fluctuation contribution to 
the thermal expansion co-efficient, ct(T), shows a clear 
SRO-induced decrease above T nem , and a sharp peak at 
T nem , as shown in Fig. [5] 

Let us now see to what extent these results square 
up with experimental data on underdoped Fe arsenides, 
where recent work on the so-called 122-systems reveals 
tantalizing signs of an orbital nematic state coupled to the 
tetragonal-orthorhombic (T-O) transition. Surprisingly, 
we find that effects of fluctuations of such an e-nematic 
order (with Ising symmetry) on the carriers, reflected in 
thermodynamic and transport data, seem to be qualita- 
tively rationalizable in our "simple" model. 

(i) the resistivity, pd c (T), shows a maximum at the T-O 
transition, not at the antiferromagnetic (AF-SDW) tran- 
sition with Q = (w, 0). It also exhibits remarkable hys- 
teretic behavior above Tt-o — 150 K. This goes hand- 
in-hand with the disappearance of the orthorhombicity, 
defined as (O) = ^r^, where a,b are unit-cell constants 
(Fe-Fe nearest neighbor distances) in the FeAs plane. In- 
terestingly, setting J' = 0, i.e, neglecting nematic or- 
dering tendency, kills the prominent peak as well as the 
hysteresis (above Tt-o) in our results, in disagreement 
with experiment (see also below). If we assume that a 
mean field nematic expectation value, (N z ) — (n/2) > 0, 
develops at the T-O transition (though it must be kept 
in mind that, because the coupling of the e-nematic or- 
der to strain, the e-nematic order cannot be separated 
from the T-O distortion, and the e-nematic transition 
will be smeared), one can address the issue of the in- 
plane resistivity anisotropy above the T-O distortion in 
our phenomenological approach. Namely, one can now 
identify the "external" zeeman field value S(T) = h z = n 
with the N z = 1/2 state with (O) > and the value 
5(T) = h z =0 with the N z = -1/2 state with (O) < in 
Ising nematic language. From our theoretical resistivity 
curves (Fig. [3]), we clearly see the development of a clear 
resistivity anisotropy: interestingly, it extends to tem- 



peratures significantly higher than Tt-o, and achieves its 
maximum precisely around Tt-o, where the nematic sus- 
ceptibility is maximal. These findings are broadly consis- 
tent with experiment^, and constitute phenomenological 
explication in terms of an electronic nematic order cou- 
pled, by symmetry, to the T-O distortion. 

To the extent that the above e-nematic-plus T-O tran- 
sition is intimately tied to a finite (O), its stabilization 
prepares the ground, via anisotropic electronic structure 
changes, for striped antiferromagnetic spin-density-wave 
(AF-SDW) state with Q = (n, 0) to emerge in a natural 
way. This program, within the context of a proposal for 
ferro-orbital order, has been carried out in Ref.@, and the 
resulting J\ a — Ju, — J2 model with large J2I Jia,b indeed 
achieves a satisfying description of spin-wave dispersion 
in the AF-SDW state in the O-structure. Electronic ne- 
maticity has hitherto not been considered within such a 
program. In view of the fact that e-nematic and T-O 
transitions are strongly coupled (the former is smeared), 
it follows that a similar instability to an AF-SDW state 
can be worked out in the present case as well. We do not 
do it here, and only mention that, once the T-O distor- 
tion occurs, the scenario of Ref. can take over. 

(ii) thermal expansion co-efficient in underdoped Fe 
arsenides has been recently measured by Wang et al£^. 
Marked anomalies in a(T) are found precisely at the 
T-O transition, and these also appear to survive and 
change with increasing doping: in particular, at doping 
levels close to the T-O boundary at low T, the "fluc- 
tuation" contribution, Aa(T), becomes negative below 
Tt-o, changing sign to positive for T > Tt-o- Re- 
markably, this is exactly the form we extract: in our 
approach (Fig. [3]), it arises due to fluctuations associated 
with Ising-like e-nematic order coupled with the T-O dis- 
tortion. Adding a purely phenomenological term, linear 
in T, to our a(T) computed above could give nice quali- 
tative agreement with these experimental results. Closer 
inspection, in fact, shows that our result for a(T) for fi- 
nite J' (i.e, including nematicity) is in much closer accord 
with data than the one with J' = 0. In particular, the 
dip (peak) in the (total) measured AVm(T)/V slightly 
below(above) Tt-o is nicely rationalized as arising from 
short-range fluctuations of an Ising-likc nematic order 
associated with (xz, yz) orbitals in real underdoped Fe 
arsenides. The estimated entropy cange across the struc- 
tural transition also bears similarities to our computed 
result, but, since the striped AF instability also occurs 
(slightly below or concurrently with) the structural in- 
stability, a direct comparison between our theory (which 
only focuses on nematicity and the accompanying T-O 
distortion) and experiment is premature. 

While the above is by no means a microscopic de- 
scription, the "localized" nematic pseudospins (N z ) can 
microscopically arise in physical situations where sizable 
multi-orbital electronic correlations selectively localize a 
subset of the relevant planar orbital band states, leav- 
ing others metallic. Thus, elucidation of the microscop- 
ies of the e-nematic-plus T-O transition in Fe-arsenides, 
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and, in particular, investigating issues like (i) how such 
e-nematic-plus T-0 transition occurs as an instability 
of the bad-metal "normal" stated—, and, (ii) its rela- 
tion (competitor) to unconventional superconductivity, 
involves much more work, and will be reported sepa- 
rately. 

Nevertheless, we close this section with a few remarks 
that may have a bearing on recent experimental data for 
underdoped Fe arsenides. Given that our phenomenolgi- 
cal Hamiltonian, 

H = t + h.c) -j'YI N * N 1 

<i,j>,o <i,j> i 

(17) 

for the e-nematic transition describes (renormalized in re- 
ality) band-like electrons coupled to a ferro- "magnetic" 
Ising model in a zeeman field, the critical behavior falls 
into the liquid-gas universality class. It is then per- 
fectly possible that an additional "tuning parameter" , 
say static chemical disorder, could tune the system to 
the T = quantum critical end-pointof the line of first- 
order transitions of this quantum liquid-gas transition. 
This could possibly bear a relation to signatures of quan- 
tum criticality observed in some of the Fe arsenides as a 
function of doping, and would not be inconsistent with 
the (observed) fact that maximal SC T c {x) in Fe ar- 
senides does not occur at the critical x = x\ where AF- 
SDW order vanishes, but rather at x — xi where Tt-o 
would have vanished, as is clear by observation of the 
T— x phase diagrams where Tt~o(x) and Tjv(x) are well- 
separated in T, x. An upshot of this reasoning would be 
then to inquire whether, in a quantum-critical scenario, 
soft fluctuations associated with an underlying QCP as- 
sociated with orbital e-nematic order could act as a pair 
glue for SC, along lines worked out by Si et a/.—. It would 
be interesting to investigate this line of thinking in more 
detail, but this is out of scope of the present work. 



V. CONCLUSION 

To conclude, in this paper, we have studied how ther- 
mally induced SRO and associated order prameter fluctu- 
ations affect various thermodynamic and resistive prop- 
erties in two cases: (i) the double exchange model with 
Jh » t, where an effective Ising- like spin model arises, 
and (ii) an electronic system undergoing a phase transi- 
tion to an e-nematic state, coupled to a lattice distortion 
and shown how a careful treatment of SRO goes quite 
a long way toward a qualitative understanding of sev- 
eral striking features in correlated systems of great cur- 
rent interest. Our analysis should be valid as an effec- 
tive phenomenological treatment in situations where the 
Ising-like order sets in independently of the nature of the 
band-like electronic state, but, once established, drasti- 
cally affects thermodynamic and transport responses via 
coupling of carriers to the order parameter spin suscepti- 
bility. Being extremely simple to implement, it can easily 
be used to analyze experimental results in a wide vari- 
ety of other systems, e.g, in multi-orbital systems where 
orbital order, generically Ising like, plays a crucial role 
in shaping their physical properties. It can also be used 
very efficiently for Ising models with competing inter- 
actions (the axial-next-nearest-neighbor Ising (ANNNI) 
model2£) in a zeeman field, which, in itself, is interesting 
as an effective model for complex ordering phenomena in 
diverse contexts. We plan to address such applications 
in future work. 
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